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Abstract 

We make use of the Nambu-Jona-Lasinio (NJL) formalism and real-time finite-temperature field 
theory to calculate hadronic current correlation functions in the deconfined phase of quantum 
chromodynamics (QCD). Here, we consider both pseudoscalar and vector currents. Since the 
method used in the lattice analysis to calculate the spectral functions requires assumptions about 
the likelihood of a particular form for the spectral function, we believe our calculations will be useful 
to researchers who wish to calculate hadronic current correlation functions at finite temperature 
using lattice-based methods. Our model makes use of temperature-dependent coupling constants 
for the NJL model. We present various arguments that such temperature dependence is necessary, if 
the results of the model are to be consistent with what is known concerning QCD thermodynamics. 
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I. INTRODUCTION 



In recent years we have seen a great deal of interest in the properties of dense matter, 
with particular attention given to diquark condensation and color superconductivity [1]. 
Since it is difficult to study the properties of dense matter in lattice simulations of QCD [2] , 
the Nambu-Jona-Lasinio (NJL) model and closely related instanton-based models have been 
used in such studies. We have become interested in a possible density dependence of the 
coupling constants of chiral Lagrangian models, since density dependence of the coupling 
parameters could affect the conclusions drawn from the studies of dense matter. We have 
introduced density-dependent coupling constants for the NJL model in earlier works [3, 
4], and have presented some arguments that such density dependence may be necessary 
[4]. However, it is much easier to discuss the temperature dependence of NJL coupling 
parameters, rather that the density dependence, since a good deal is known concerning 
finite-temperature QCD thermodynamics [5]. In particular, the study of a gluon gas at high 
temperature suggests that the system is well described as weakly interacting for T > 6T C , 
where T c is the temperature of the confinement- deconfinement phase transition. (See Fig. 1.3 
ofRef. [5].) 

In this work we suggest that rather straightforward calculations of hadronic current 
correlation functions at finite temperature can provide information concerning a possible 
temperature dependence of the NJL coupling parameters. We perform calculations of 
hadronic current correlation functions for pseudoscalar and vector currents in the range 
1.2 < T/T c < 5.88. We make use of two models for the temperature dependence of the 
NJL coupling parameters. For model 1, we use G(T) = G[l — 0.17(T/T C )], which was the 
form used in our previous studies of meson properties [6] and hadronic current correlation 
functions [7] at finite temperature. In this work we also introduce a model 2, for which 
G(T) = G[l - 0.0289(T/T C ) 2 ]. In both cases, G(T) = for T/T c = 5.88. For values of 
T/T c > 5.88 we put G(T) = for both models. 

For the sake of completeness we present the Lagrangian of a generalized NJL model that 
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we have used in our studies of meson properties at finite temperature and density 
C = q(i$> - m°)q + ^ £[(gA\z) 2 + (qi^Xq) 2 ] 

i=0 

-^X>aW + (^w) 2 ] 

i=0 

+ ^{det[g(l + 75 )g] + det[g(l - 75 )g]} 

+C con f . (1.1) 

Here, m° is a current quark mass matrix, m° = diag(m°, ra°, m°). The A* are the Gell- 
Mann (flavor) matrices. Here, A = a/2/31 with 1 being the unit matrix. The fourth term 
on the right-hand side of Eq. (1.1) is the 't Hooft interaction. Finally, £ co „/ represents the 
model of confinement we have used in our work. 

In order to specify the parameters in Eq. (1.1) we may refer to Ref [8]. There, a quite 
detailed fit was made to the properties of the 77 mesons. Our analysis yields fits the first 
four experimentally known levels and provides excellent fits to the mixing angles and decay 
constants of the 7/(547) and 77/(958). That work led to the specification of G,s = 11.84 
GeV -2 . The 't Hooft interaction strength, Gd, was taken to be in the range —220 GeV~ 5 
< Gd < —180 GeV -5 when calculating the properties of the 77 mesons. We also used 
Gy = 13.0 GeV -2 in Ref. [8]. (It is worth noting that Hatsuda and Kunihiro [9] have used a 
sharp moment cutoff of A 3 = 0.6314 GeV and a value of Gd = 185.1 GeV -5 . Their value of 
Gs is smaller than ours, since they do not include a model of confinement in their analysis.) 

In terms of these parameters we can introduce effective coupling constants for states 
with particular quantum numbers. The effective coupling parameter involves the vacuum 
condensates a =< Q\uu\Q >, (3 =< 0\dd\0 > and 7 =< 0\ss\0 >. For the pseudoscalar- 
isovector states, we have [9] 

Gp ff = G s + ^. (1.2) 

We take a = (3 = -(0.241GeV) 3 and 7 = -(0.258GeV) 3 [10]. If we put G D = -190 GeV~ 5 , 
we find G P = 13.47 GeV~ 2 . 

In Ref. [6] we presented relativistic random phase approximation (RPA) calculations for 
the pion, kaon and several other mesons. It was found that in the case of pion, that the use 
of G e P ff = 13.49 GeV" 2 gave a pion energy of 140 MeV at T = 0. It is gratifying to see that 
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an independent calculation gives rise to essentially the same effective coupling constant as 
that obtained from the parameters determined in Ref.[8]. 

If we consider the scalar-isovector states, the effective coupling constant is 

G e / f = G s -^, (1.3) 

which, for G s = -11.84 GeV -2 , G D = 190.0 GeV -5 and 7 = -(0.258GeV) 3 , yields G s = 
10.21 GeV -2 . When we compare G e / f with G e /* , we see that we can expect the resonant 
enhancement of the spectral functions for the pseudoscalar-isovector channel will be larger 
in magnitude than that of the spectral function of the scalar-isovector channel. On the 
other hand, since the 't Hooft interaction does not affect the vector or axial-vector coupling 
constants, these coupling constants would be equal in our model, leading us to expect quite 
similar spectral functions for these two channels. 

In our study of the 77 mesons we used Gy = 13.00 GeV -2 . However,the results for the 77 
mesons are not particularly sensitive to that parameter. In earlier work, with a sharp cutoff 
for the NJL model of A 3 = 0.622 GeV, we used G v = 12.46 GeV -2 to fit the energy of the 
o;(782) [11]. In other calculations we have used Gy = 11.46 GeV -2 . In this work we present 
the results for the latter value. (We have checked that the spectral functions calculated with 
either Gy = 12.46 GeV -2 or G v = 11.46 GeV -2 differ by less than 4 percent.) 

It was found in our various studies of meson properties, made after we completed Ref. 
[11], that the use of a sharp cutoff did not allow us to study radial excitations with a 
large number of nodes. However, we found that a Gaussian regulator, exp[— k 2 /a 2 }, with 
a = 0.605 GeV, gave similar results as that of the sharp cutoff A 3 = 0.622 GeV for the low- 
lying states and also allowed us to describe radial excitations with many nodes. Therefore, 
all our more recent calculations have been made using the Gaussian regulator with a = 0.605 
GeV. (The coupling parameters obtained in Ref. [8], for example, were determined with this 
Gaussian regulator, so that any change of the regulator parameter, a, would require that 
we provide a new set of coupling parameters corresponding to the new value of a.) 

We have recently reported results of our calculations of the temperature dependence of 
the spectra of various meons [6]. These calculations were made using our generalized NJL 
model which includes a covariant model of confinement. We have presented results for the 
71, K, do, fo and Kq mesons in Ref. [6]. In that work, temperature-dependent constituent 
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quark masses were calculated using the equation [12] 

m(T) = m° + 2G S (T)N C ^P- I dp^ tanh(lf^) . (1.4) 

Here, m° is the current quark mass, G$(T) is a temperature-dependent coupling constant 
introduced in our model. (In Ref.[6] we used model 1 for Gs(T).) Here, N c = 3 is the 
number of colors, f3 = 1/T and E p = [p 2 + m 2 (T)] 1/2 . Further, A = 0.631 GeV is a cutoff 
such that \p\ < A. We have put A = 0.631 GeV, since that is the cutoff that is often used 
when solving Eq. (1.4) [9]. We have also used mJJ = 5.5 MeV and m° s = 120 MeV. Since in 
our work dealing with meson spectra, we have used m u = 0.364 GeV and m s = 0.565 GeV 
as phenomenological parameters, we have put G$(0) = 2G, with G = 5.691 GeV -2 , in the 
notation of Ref.[12], to reproduce those values. In Ref.[12] G is one-half of the Gs defined 
in Eq.(l.l). 

Since the use of temperature-dependent coupling constants is an unusual feature of our 
work, we now describe some advantages that follow from that choice. Our original intro- 
duction of temperature dependence, G(T) = G[l — 0.17(T/T C )], with T c = 170 MeV, had 
the purpose of simulating the dynamical interactions which eliminate pion or kaon conden- 
sation at relatively low temperatures. However, we found when solving Eq.(1.4) for the 
temperature-dependent constituent mass of the up quark, chiral symmetry was (partially) 
restored at lower temperatures than when constant coupling parameters were used. This 
may be seen in Fig.l, where we show the results that follow from the use of a constant 
value of G(T) = 5.691 GeV~ 2 (in the notation of Ref.[12]) as a dotted line. Here we have 
used m° = 5.5 MeV as the current quark mass. For the solid and dashed curves in Fig.l 
we have used G(T) = G[l — 0.17(T/T C )]. For the sake of this discussion, let us consider 
the reduction of the constituent mass to 50 MeV from 364 MeV as signal of the (partial) 
restoration of chiral symmetry. For the dashed and solid curves, that restoration take place 
at about T c = 170 MeV, while for the dotted curve the (partial) restoration of chiral sym- 
metry takes place at about 250 MeV. For QCD, the transition temperature is about 150-170 
MeV. Since we are attempting to create a model that has some correspondence to QCD, 
with dynamical quarks, we see that the temperature-dependent coupling constants lead to 
the desired behavior. 

The temperature dependence of the constituent mass in the case of temperature- 
dependent coupling constants is such that we described the mesonic confinement- 
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deconfinement transition as taking place at T « T c with T c = 170 MeV in Ref. [6]. In that 
work we have studied the confinement- deconfinement transition for the n, K, ao, fo and Kq 
mesons. We have checked that, without the introduction of the temperature-dependent cou- 
pling constants, the confinement-deconfinement transition would take place at a significantly 
higher temperature. 

A further advantage of the use of our temperature-dependent coupling parameter is that 
the NJL interaction goes to zero for T = 5.88T c .(It is put equal to zero for T > 5.88T C .) We 
suggest that that is consistent with QCD thermodynamics, since it appears that the system 
is weakly coupled for T > 6T C [5]. 

In Fig. 2 we show the behavior of both the up and strange quark constituent masses 
calculated with Eq.(1.4) using current masses mP u = 5.50 MeV and m Q s = 120 MeV. For 
the calculations made in this work, we have used the temperature-dependent up (or down) 
quark masses shown in Fig.2 for temperatures for which the quark masses are quite small. 
Therefore, our results are insensitive to variations of the mass values shown in Fig.2. 

In calculating the constituent mass values we have neglected the confining interaction. 
That interaction was taken into account in our earlier Euclidean-space calculation of the 
quark self-energy [13], which also included the effects related to the 't Hooft interaction. We 
found that, to a good approximation, we could neglect the confining and 't Hooft interactions, 
if we modified the value of the NJL coupling constant, Gs, and we adopt that approach when 
using Eq.(1.4). 

The organization of our work is as follows. In Sec. II we review our calculations of po- 
larization functions at finite temperature. In Sec. Ill we discuss the calculation of hadronic 
current correlation functions, making use of the results presented in Sec. II. In Sections IV 
and V we present the results of our numerical calculations of correlators of pseudoscalar and 
vector currents, respectively. We also compare our results to some recent lattice calculations 
of pseudoscalar and vector correlators [14-16]. Finally, in Sec. VI we present same further 
discussion and conclusions. 

II. POLARIZATION FUNCTIONS AT FINITE TEMPERATURE 

The basic polarization function that is calculated in the NJL model is shown in Fig. 3. We 
will consider calculations of such functions in the frame where P = 0. In our earlier work, 
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FIG. 1: We exhibit the solution for m u (T) obtained using Eq.(5.38) of Ref.[12], where m° = 5.50 
MeV and A = 0.631 GeV. The dotted curve corresponds to the use of a constant value G = 5.691 
GeV -2 , in the notation of Ref. [12]. For the solid and dashed curves we have used G(T) = 
G[l — 0.17(T/T C )]. For the solid curve T c = 0.170 GeV and for the dashed curve we have used 
T c = 0.150 GeV. 

calculations were made after a confinement vertex was included. That vertex is represented 
by the filled triangular region in Fig. 3. However, we here consider calculations for T > 1.2 T c 
where confinement may be neglected. We will, however, use the temperature-dependent mass 
values shown in Fig. 2. 

The procedure we adopt is based upon the real-time finite-temperature formalism, in 
which the imaginary part of the polarization function may be calculated. Then, the real 
part of the function is obtained using a dispersion relation. The result we need for this work 
has been already given in the work of Kobes and SemenofT [17]. (In Ref. [17] the quark 
momentum in Fig.3 is k and the antiquark momentum is k — P. We will adopt that notation 
in this section for ease of reference to the results presented in Ref. [17].) With reference to 
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FIG. 2: Temperature dependent constituent mass values, m u (T) and m s (T), calculated in a mean- 
field approximation [12], are shown. [See Eq. (1.4)]. Here m° u = 0.0055 GeV, m° s = 0.120 GeV, 
and G(T) = 5.691 [1 - 0.17(T/T C )], if we use Klevansky's notation [12]. (The value of G s of Eq. 
(1.1) is twice the value of G used in [12].) 




FIG. 3: The upper figure represents the basic polarization diagram of the NJL model in which 
the lines represent a constituent quark and a constituent antiquark. The lower figure shows a 
confinement vertex [filled triangular region] used in our earlier work. For the present work we 
neglect confinement for T > 1.2 T c , with T c = 150 MeV. 
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Eq. (5.4) of Ref. [17], we write the imaginary part of the scalar polarization function as 

W S (P>,X) = l m)MP ^ e -^ (2.D 

{(1 - m(A;) - n 2 (k))5(P° - E x {k) - E 2 (k)) 
-(ni(k) - n 2 (k))5(P° + E 1 (k) - E 2 {k)) 
-(n 2 (k) - ni(k))5(P° - E ± (k) + E 2 {k)) 
_(1 _ ni (k) - n 2 (k))5(P° + E 1 (k) + E 2 (k))} . 

Here, Ei(k) = [k 2 + m^T)] 1 / 2 . Relative to Eq. (5.4) of Ref. [17], we have changed the 
sign, removed a factor of g 2 and have included a statistical factor of 2N C , where the fac- 
tor of 2 arises from the flavor trace. In addition, we have included a Gaussian regulator, 
exp[— k 2 /a 2 ], with a = 0.605 GeV, which is the same as that used in most of our applications 
of the NJL model in the calculation of meson properties. We also note that 

"*(*) = jnm + 1 » W 

and 

n2 ^ = + x • ( 2 - 3 ) 

For the calculation of the imaginary part of the polarization function, we may put k 2 = 
m\(T) and (k — P) 2 = m|(T), since in that calculation the quark and antiquark are on- 
mass-shell. In Eq. (2.1) the factor (3 S arises from a trace involving Dirac matrices, such 
that 

(3 S = -TrP + rriiXJf - f + m 2 )\ (2.4) 
= 2P 2 -2(m 1 + m 2 ) 2 , (2.5) 

— * 

where mi and m 2 depend upon temperature. In the frame where P = 0, and in the case 
mi = m 2 , we have (5s = 2P 2 (1 — 4m 2 / Pq). For the scalar case, with mi = m 2 , we find 

Im J S {P\ T) = ^{l- A -^) ^ e- 2 /* 2 [1 - 2^*)] , (2.6) 

where 

k 2 = ^-m 2 (T). (2.7) 
9 



For pseudoscalar mesons, we replace f3s by 

P P = -Tr[i l5 (ft + m 1 )i l5 (lk-P + m 2 )} (2.8) 
= 2P 2 - 2(mi - m 2 f , (2.9) 

which for m x = m 2 is (3 P = 2P 2 in the frame where P — 0. We find, for the 7r mesons, 

Im Jp(P 2 ,T) = ^ (l - ^)' /2 e- J V-'[l - 2» lW ] , (2.10) 

where k 2 = P 2 /4 — m 2 (T), as above. Thus, we see that, relative to the scalar case, the 
phase space factor has an exponent of 1/2 corresponding to a s-wave amplitude. For the 
scalars, the exponent of the phase-space factor is 3/2, as seen in Eq. (2.6). 

For a study of the vector-isovector correlators, we introduce conserved vector currents 
3iiA x ) = <?( x )o 'n^iQ( x ) with i=l, 2 and 3. In this case we define 



and 



jp(P^T)=[g^-^^)j v (P 2 ,T) 



C^(P 2 ,T) =[g^- ) C V (P 2 ,T) , 



(2.11) 



(2.12) 



taking into account the fact that the current j^i(x) is conserved. We may then use the fact 
that 



Jv(P 2 ,T) = - 9iiu JP(P 2 ,T) 



(2.13) 



and 



lmJ v (P 2 ,T) 



P 2 + 2m 2 u (T) 
An 



-lmJ P (P 2 ,T) 



(l-^) 1/2 e^ 2 [l-2n lW ] (2.14) 

(2.15) 



We now consider 



C = Tr[ 7 ^ + mi) 7l/ (# - P + m 2 )\ , 



(2.16) 



and calculate 



<T0L = 4[P 2 - m\ - m\ + 4m im2 ] , 



(2.17) 
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which, in the equal- mass case, is equal to 4P 2 + 8m 2 (T), when m x = m 2 and P — 0. This 
result will be needed when we calculate the correlator of vector currents in the next section. 
Note that for the elevated temperatures considered in this work m u (T) = m d (T) is quite 
small, so that 4Pq + 8m 2 (T) can be approximated by 4Pq when we consider the vector 
current correlation functions. In that case, we have 

1m J V (P 2 , T) ~ hm J P (P 2 ,T), (2.18) 

At this point it is useful to define functions that do not contain the Gaussian regulator: 

Im J P (P 2 ,T) = ^ (l " A -^Y [1 - 2n!(*)] , (2.19) 

and 

Im J V (P 2 , T) = H^JL (l - V2 [1 - 2n 1 (Jb)] . (2.20) 

For the functions defined in Eq. (2.19) and (2.20) we need to use a twice-subtracted disper- 
sion relation to obtain Re Jp(P 2 , T), or Re Jy(P 2 , T). For example, 

Re J P (P 2 , T) = Re J P (0, T) + [Re J P (P 2 , T) - Re J P (0, T)\ + (2.21) 

M) 

P 2 (P 2 ~ Pj) [ R2 Im J P ( S ,T) 

L 2( T) %(p 2 - S )(p 2 - S )' 

where A 2 can be quite large since the integral over the imaginary part of the polarization 
function is now convergent. We may introduce Jp(P 2 ,T) and Jy(P 2 ,T) as complex func- 
tions, since we now have both the real and imaginary parts of these functions. We note that 
the construction of either Re Jp(P 2 ,T) or Re Jy(P 2 ,T) by means of a dispersion relation 
does not require a subtraction. We use these functions to define the complex functions 
J P (P 2 ,T) and Jy(P 2 ,T). 

In order to make use of Eq.(2.21) we need to specify ReJp(O) and ReJp(P 2 ). We found it 
useful to take P 2 = -1.0 GeV 2 and to put ReJ P (0) = ReJ P (0) and ReJ P (P 2 ) = ReJ P (P 2 ). 
The quantities ReJy(O) and ReJy(PQ) are determined in an analogous function. This proce- 
dure in which we fix the behavior of a function such as J P (P 2 ) or Jy(P 2 ), which may be used 
when making calculations for large P 2 , is quite analogous to the procedure used in Ref. [18]. 
In that work we made use of dispersion relations to construct a continuous vector-isovector 
current correlation function that had the correct perturbative behavior for P 2 — > — oo and 
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also described the low-energy resonance present in the correlator due to the excitation of the 
p meson. In Ref. [18] the NJL model was shown to provide a quite satisfactory description 
of the low-energy resonant behavior of the vector-isovector correlation function. 

We may compare our expressions for (l/7r)ImJ(P 2 ) with those given in Eq.(2.7) of 
Ref. [19]. We find that the scalar and pseudoscalar polarization functions are defined there 
with a sign opposite to ours. However, since we have defined the pseudoscalar current 
using the matrix 275 rather than 75, which was used in Ref. [19], our expression for the 
spectral function agrees with that of Ref. [19] with respect to the sign in the pseudoscalar 
case. In the case of the vector and axial-vector spectral functions, we have the same sign 
convention as that of Ref. [19]. We disagree with the expressions for the spectral functions 
given for the scalar and axial- vector currents in Ref. [19] with respect to the exponent of the 
phase-space factor. In Ref. [19], the phase-space factor for the scalar and axial- vector case is 
[1 — (2m/w) 2 ]^ 2 , while we have [1 — (2m/cj) 2 ] 3 / 2 for those spectral functions. We agree that 
the phase-space factor is [1 — (2m/w) 2 ]^ 2 for the pseudoscalar and vector current spectral 



In Ref. [19] use is made of the hard-thermal-loop (HTL) approximation. It was found that 
the scalar channels are only moderately influenced by the HTL medium effects, while the 
HTL vertex corrections lead to divergent vector correlators. 

III. CALCULATION OF HADRONIC CURRENT CORRELATION FUNCTIONS 

In this section we consider the calculation of temperature-dependent hadronic current 
correlation functions. The general form of the correlator is a transform of a time-ordered 
product of currents, 



where the double bracket is a reminder that we are considering the finite-temperature case. 

For the study of pseudoscalar states, we may consider currents of the form jp,i(x) = 
q(x)i , y^X 1 'q(x), where, in the case of the 7r mesons, % — 1,2 and 3. For the study of scalar- 
isoscalar mesons, we introduce js,i( x ) — q{x)\ l q(x), with % = for the flavor-singlet current 
and i = 8 for the flavor-octet current [7]. 

In the case of the pseudoscalar-isovector mesons, the correlator may be expressed in terms 



functions. 




(3.1) 



12 



of the basic vacuum polarization function of the NJL model, Jp(P 2 ,T) [9,10,12]. Thus, 

1 



C P (P\T) = Jp(P\T) — 



(3.2) 



G P {T)Jp{P\T)' 

where Gp(T) is the coupling constant appropriate for our study of n mesons. We have 
found G P (T) = 13.49 GeV~ 2 by fitting the pion mass in a calculation made at T = 0, with 
m u = m d = 0.364 GeV. The result given in Eq. (3.2) is only expected to be useful for small 
P 2 , since the Gaussian regulator strongly modifies the large P 2 behavior. Therefore, we 
suggest that the following form is useful, if we are to consider the larger values of P 2 



C P (P 2 ,T) 
P 2 



Jp{P\T) 
P 2 



1 



l-G P (T)Jp(P 2 ,T) 



(3.3) 



(As usual, we put P = 0.) This form has two important features. At large Pq, 
lmCp(P ,T)/P^ is a constant, since Im J p (Pq,T) is proportional to Pq. Further, the de- 
nominator of the second term on the right-hand side of Eq.(3.3) goes to 1 for large Pq. On 
the other hand, at small Pq, that denominator is capable of describing resonant enhance- 
ment of the correlation function. (We may again refer to Ref. [18], in which a similar 
approximation is described.) 
We then have 



C V {P 2 ,T) = J V (P 2 ,T)- 



1 



G V (T)J V (P*,T)' 



(3.4) 



where we have introduced 



ImJ v (P 2 ,T) = - 



P 2 + 2m 2 u (T) 



47T 



^P) 1/2 [l-2„ lW1 



~ -lmJ P (P 2 ,T). 



(3.5) 
(3.6) 



(See Eq.(2.7) for the specification of k — \k\.) In the literature, u is used instead of Pq [14 
-16]. We may define the spectral functions 



1 

cry (u, T) = — Im C v (co, T) , 

71 



(3.7) 



and 



a P (uj,T) = -hnC P (u,T). 

TT 



(3.8) 
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We may use the notation ap{uj,T) and ay(u;,T) for the spectral functions given in the 
literature [14-16]. We have the following relations: 



a P (co,T) = a P (u,T) (3.9) 

and 

a v (u,T) 3 

2 = ^v{u,T), (3.10) 

where the factor of 3/4 arises because, in Refs. [14-16], there is a division by 4, while we 
have divided by 3 as in Eq. (2.13). 

In addition to providing values of oy(cj, T)/u 2 and <Jp(u, T)/u 2 , the r-dependent 
hadronic correlators Gp{r,T) and Gy{r,T) were introduced [14-16], with 



G(t,T)= / daJcr(cJ,T)fr(cc>,r,T). (3.11) 
Jo 



Here, 



cosh[c<j(r — ] 

K{ ^= kkg) - (3 ' 12) 

We will present results of our calculations of G P (r,T) and G v (t,T) in Sections IV 
and V, respectively. We will also present the values obtained for Im Cp(P , T)/P 2 and 
ImCy(Po) T)/Pq for various values of T/T c in Sections IV and V. 

We note that values have been obtained for the spectral functions in the scalar, vector, 
pseudoscalar and axial- vector channels in Ref . [20] . These functions exhibit resonant features 
at about 2, 7 and 14 GeV when T ps IAT C . High-energy resonances are still present at 5 and 
13 GeV, when T 1.9T c . The interpretation of these high-energy resonances is uncertain 
[21]. 



IV. CORRELATION FUNCTIONS FOR PSEUDOSCALAR-ISOVECTOR CUR- 
RENTS - NUMERICAL RESULTS 

In Fig. 4 we show lmC P (P ,T)/ P 2 for T/T c = 1.5, and for model 1, as a solid line. The 
dashed line represents the result at T/T c = 1.5, if we put G P (T) = 0. The results for a 
broader range of temperatures are shown in Fig. 5. In Fig.6 we show the ratios Rp(P , T) = 
lmC P {P ,T) /lmC { p ] {P ,T) calculated for various temperatures. Here, lmCP{P ,T)/P 2 is 

14 



5 10 15 20 25 

P „ /T 

c 

FIG. 4: Values of ImC P (P , T)/P$ are shown at T/T c = 1.5 for model 1, where G P (T) = G P [1 - 
0.17 (T/T c )] (solid line). The dashed line represents the result obtained when T/T c = 1.5 and 
G P {T) = 0. (Here, the dashed line represents the values of Im J P (P , T)/P 2 for T/T c = 1.5.) 

calculated with G P (T) = 0. In Fig.7 we show the ratio R ( a\r,T) = G P (r,T) /G { p\r,T) 
for various temperatures. Here, Gp\r,T) represents G p (t,T) evaluated for G P (T) = 0. In 
Fig. 8 we exhibit Im C P (P , T)/P 2 for various temperature for model 2. 

In Fig. 9 we show, as a dashed line, values of Im C P (P , T)/P 2 at T/T c = 5.88 obtained 
using a constant value of G P (T) = 13.49 GeV -2 . The solid line is the result for model 1 (or 
model 2) at T/T c = 5.88, in which case G P (T) = for both models. Since we have argued 
that for T/T c ~ 6 the quark-gluon plasma should be a weakly interacting system, the use of 
a constant value for G P (T) appears to be unacceptable. We present this result to support 
our argument that the coupling parameters of the chiral Lagrangian model should be made 
temperature dependent to be consistent with QCD thermodynamics. 

In Fig. 10 we show values of G p (t,T)/T 3 for T/T c = 1.5 [dashed line], T/T c = 3.0 [dash- 
dotted line] and T/T c = 5.88 [solid line]. In Fig.l of Ref. [11] we see a marked difference 
in the behavior of G p (t,T)/T 3 and Gy(r, T)/T 3 . (We will present some discussion of this 
matter in Sec.V.) 
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FIG. 5: Values of Im Cp(Pq,T)/Pq are shown for model 1 and for various temperatures: T/T c = 1.2 
[sold line], T/T c = 1.5 [dashed line], T/T c = 2.0 [dotted line], T/T c = 3.0 [dashed-dotted line], 
T/T c = 4.0 [dashed-double dotted line], and T/T c = 5.88 [short-dashed line]. 




FIG. 6: The ratio R P (P ,T) = lmC P (P ,T) /ImC},' (P ,T) is shown for various values of T/T c = 
1.2 [solid line], T/T c = 1.5 [dashed line], T/T c = 2.9 [dotted line], T/T c = 3.0 [dot-dashed line] 
and T/T c = 4.0 [double dot-dashed line]. Here, lmCp\Po,T) is calculated at temperature T with 
G P {T) = 0. 
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FIG. 7: Values of R { ^'{t,T) are shown for values of T/T c = 1.2 [solid line], T/T c = 1.5 [dashed 
line], T/T c = 2.0 [dotted line], T/T c = 3.0 [dashed-dotted line], and T/T c = 4.0 [dashed-double 
dotted line]. For T/T c = 5.88, R { p(T,T) = 1. 

V. CORRELATION FUNCTIONS FOR VECTOR-ISOVECTOR CURRENTS - 
NUMERICAL RESULTS 

In Fig. 11 we show the values of Im C v {Po, T)/Pq for model 1 and for T/T c = 1.5, as 
a solid line. The result for T/T c = 1.5 and GV(T) = is represented by the dashed line. 
[See Fig. 4.] In Fig. 12 we show values of Im Cy(P , T)/Pq for various values of T/T c and for 
model 1. Corresponding results for model 2 are given in Fig. 13. 

In Fig.14 we show the ratio R V (P ,T) = lmC v {P ,T)/lmC { y\Po,T), where 
lmC$\P ,T) is calculated with G V (T) = 0. (See Fig.9). In Fig.15 we show the ratio 
R^\t,T) = G v (T,T)/G { y\r,T), where G v (0)(r,T) represents G v (r,T) calculated with 
G V (T) = 0. In Fig.16 we show, as a dashed line, values for T/T c = 5.88 of lmC v {P , T)/P 2 , 
which were calculated with a constant value of Gy(T) = 11.46 GeV~ 2 . The solid line repre- 
sents the results for G V {T) = 0, which is characteristic of models 1 and 2 when T/T c = 5.88. 
The comments made with respect to Fig.9 are also applicable here. 

In Fig.17 we show values of G v {r,T)/T 3 for T/T c = 1.5 [dashed line], T/T c = 3.0 [dash- 
dotted line] and T/T c = 5.88 [solid line]. These results were obtained with model 1, and 
may be compared to those shown in Fig.9, where we see generally similar behavior. That 
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FIG. 8: Values of lmC P (Po, T)/Pq are shown at for various temperatures for model 2, with 
G P (T) = G P [1 - 0.0289 (T/T c ) 2 ]. (See the caption to Fig.5.) 




FIG. 9: Values of Im C P (P , T)/P 2 are shown for T/T c = 5.88. The dashed line represents the 
result when G P {T) = G P = 13.49 GeV" 2 . The solid line is the result for G P {T) = 0, which is 
characteristic of models 1 and 2, when T/T c = 5.88. (The solid line, therefore, represents the values 
of Im J P (P ,T)/P$ for T/T c = 5.88.) 
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FIG. 10: Values of G p (t, T)/T 3 are shown for T/T c = 1.5 [dashed line], T/T c = 3.0 [dashed-dotted 
line], and T/T c = 5.88 [solid line]. These results were calculated with model 1. 

0.20 -, 




FIG. 11: Values of ImC v {P ,T)/P$ are shown for T/T c = 1.5. (See the caption of Fig.4.) 

is in strong contrast to the results shown in Ref. [15], where Gy{r,T) /T 3 and Gp(r,T)/T 3 
shows quite different behavior, with the result for the vector correlator close to the values 
for G(T) = at T/T c = 1.5 and T/T c = 3.0. This suggests that the value of G v = 11.46 
GeV -2 that we have used in this work may be too large, or that the temperature dependence 
of Gy(T) is such as to yield smaller values than those obtained in this work for model 1 
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FIG. 12: Values of ImCv{Po,T) / Pq are shown for model 1 and for various temperatures. (See the 
caption of Fig. 5.) 

0.30- 




FIG. 13: Values oilmCy(Po,T) / Pq are shown for model 2 and for various temperatures. (See the 
caption of Fig. 5.) 

or model 2. (In general, we prefer the the second of these possibilities, since our value 
of Gy at T = is set by fitting meson spectra.) In Fig. 18 we show Im Cy(P , T)/Pq for 
T = 1.5 T c , for model 1, with Gy = 11.46 GeV~ 2 [solid line], G v = 8.00 GeV" 2 [dashed 
line], G v = 4.0 GeV~ 2 [dotted line] and Gy = 0.0 GeV' 2 [dot-dashed line]. It would be of 
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FIG. 15: Values of Rq (t, T) are shown for values of T/T c = 1.2 [solid line], T/T c = 1.5 [dashed 
line], T/T c = 2.0 [dotted line], T/T c = 3.0 [dashed-dotted line], and T/T c = 4.0 [double dot-dashed 
line]. For T/T c = 5.88, R ( p(r,T) = 1. 

interest to consider the results of Ref. [15] for the Euclidean-space correlator to be correct 
and to then determine GV(T) so that we fit the values of GV(r, T)/T 3 obtained in the lattice 
calculations. We defer such a project to a future work. 



21 



5 10 15 20 25 

P „ /T 

c 

FIG. 16: Values of Im C V (P , T)/P 2 are shown for T/T c = 5.88. The dashed line represents the 
result when G V (T) = Gy = 11.46 GeV~ 2 . The solid line is the result for G V (T) = 0, which is 
characteristic of models 1 and 2, when T/T c = 5.88. The solid line, therefore, represents the values 
of Im J V (P ,T)/P$ for T/T c = 5.88. 
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FIG. 17: Values of G v (t, T)/T 3 are shown for T/T c = 1.5 [dashed line], T/T c = 3.0 [dashed-dotted 
line], and T/T c = 5.88 [solid line]. These results were calculated with model 1. 
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FIG. 18: Values of ImC v (P , T)/P$ are shown for T/T c = 1.5 with G v = 11.46 GeV" 2 [solid line], 
G v = 8.00 GeV~ 2 [dashed line], G v = 4.00 GeV" 2 [dotted line] and G v = 0.0 GeV~ 2 [dot-dashed 
line]. Here, we have used model 1. 

VI. DISCUSSION 

It is difficult to make a definitive comparison of our results and the results obtained 
for the spectral functions in the lattice simulations, since those results are accompanied by 
large errors. Also, it is somewhat difficult to understand why the different behavior seen for 
G v (t,T)/T 3 and G p (t,T)/T 3 in Ref. [15], leads to values for ct p {uj,T)/lu 3 and a v {u,T)/u 3 
that are rather similar [14-16]. On the other hand, the chiral Lagrangian model provides a 
systematic study, which may yield some guidance for further studies of lattice QCD. This 
may be particularly important in the light of the comments made in Ref. [15]: 

"The reconstruction of a H {u,T) and in particular the determination of its low 
energy structure thus is difficult at non-zero temperature. Additional compli- 
cations arise in lattice calculations which necessarily are performed on lattices 
with finite number of points (N T ) in Euclidean time. The correlation functions 
Gh(t, T) can thus be calculated only at finite set of Euclidean times tT = k/N T , 
with k — 0, ...N T — 1. In order to reconstruct the spectral functions from this 
limited set of information it is necessary to include in the statistical analysis 
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of numerical results also prior information on the structure of Gh{t, T) as well 
as assumptions about the likelihood of a certain spectral function an{^,T). It 
has been suggested to provide this additional information through the applica- 
tion of the Maximum Entropy Method (MEM) [22, 23], which has been applied 
successfully to many other ill-conditioned problems in physics • • • ." 

While it may be premature to compare our results for the spectral functions with those 
given in literature, the errors for the Euclidean-time correlation functions are small. There- 
fore, in Ref. [24] we have made an extensive comparison of our results with those presented 
in Ref. [15]. In that work we provided additional evidence for the temperature dependence 
of the coupling parameters that we have introduced. 
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